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ABSTRACT 


We report high-precision multi-wavelength linear-polarization observations of the bright B9 (or A0) 
star e Sagittarii. The polarization shows the distinctive wavelength dependence expected for a rapidly 
rotating star. Analysis of the polarization data reveals an angular rotation rate w (= Q/Qerit) of 0.995 
or greater, the highest yet measured for a star in our galaxy. An additional wavelength-independent 
polarization component is attributed to electron scattering in a low-density edge-on gas disk that 
also produces the narrow absorption components seen in the spectrum. Several properties of the star 
(polarization due to a disk, occasional weak Ha emission, and multiple periodicities seen in space 
photometry) resemble those of Be stars, but the level of activity in all cases is much lower than that of 
typical Be stars. The stellar properties are inconsistent with single rotating-star evolutionary tracks, 
indicating that it is most likely a product of binary interaction. The star is an excellent candidate 
for observation by interferometry, optical spectropolarimetry to detect the Óhman effect, and UV 
polarimetry; any of which would allow its extreme rotation to be tested and its stellar properties to 


be refined. 
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1. INTRODUCTION 


'The non-spherical shape of a rapidly rotating star re- 
sults in a net linear polarization of the integrated light 
of the star. This effect was first predicted by Harring- 
ton & Collins (1968), but the resulting polarizations 
are small, and have only recently been observed (Cot- 
ton et al. 2017a; Bailey et al. 2020a; Lewis et al. 2022; 
Howarth et al. 2023). The polarization has a distinctive 
wavelength dependence, and is strongly dependent on 
the rotation rate w, the axial inclination and the grav- 
ity of the star. This allows polarization data, combined 
with other observations, to set significant constraints on 
the stellar properties. 
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€ Sagittarii (HD 169022, hereafter e Sgr) is a bright 
(V = 1.85) rapidly rotating star at a distance of 44 pc 
(van Leeuwen 2007) and with a number of anomalous 
properties. The spectral type is variously reported as 
B9IVp (Slettebak 1975), B9.5III (Houk 1982), AOI 
(Gray & Garrison 1987) and AOIII (Paunzen et al. 
2001). Narrow, "shell" absorption components in the 
Balmer and other lines are seen (Slettebak 1975), and 
Gray & Garrison (1987) call it a “proto-shell star". It 
has been suggested to be a A Boo star (Slettebak 1975). 
However, Baschek & Slettebak (1988) conclude, from a 
detailed abundance study, that e Sgr, while “peculiar 
in both spectrum and colors, is definitely not a A Boo 
star". 

A companion star of magnitude V ~ 7.7 at a dis- 
tance of ~ 2.1 arc seconds was discovered by Golimowski 
et al. (1993). Infrared observations by Hubrig et al. 
(2001) give the K magnitude of the companion as 6.5, 
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and are consistent with it being a main-sequence star of 
0.95 Mo with Teg= 5808 K. 

€ Sgr was found to have an infrared (IR) excess in 
60 um Infra Red Astronomical Satellite (IRAS) data 
by Cote (1987). The IRAS excess was interpreted by 
Rhee et al. (2007) as originating from a dust disk with 
a temperature of 100 K at a radius of 3.5 arc seconds 
from the star (equivalent to 150 au) and a fractional 
luminosity (Laust/L. = 4.46 x 1079). The IR excess 
was also detected by Spitzer at 13 and 31 um (Chen 
et al. 2014). 

Polarization observations of e Sgr were reported by 
Cotton et al. (2016a) and Bailey et al. (2017), and were 
larger than expected at its distance, implying the pres- 
ence of intrinsic polarization. In this paper we report 
a more extensive set of polarization observations and 
compare these with models of rapidly rotating stars. 


2. OBSERVATIONS 
2.1. Polarization Measurements 


We report here 28 linear polarization observations of 
€ Sgr. Of these, 26 were made on the 3.9-m Anglo- 
Australian Telescope (AAT) at Siding Spring Observa- 
tory using either the HIPPI (high precision polarimetric 
instrument, Bailey et al. 2015) or HIPPI-2 (Bailey et al. 
2020b) polarimeters. Two observations were made with 
the Mini-HIPPI (Bailey et al. 2017) polarimeter on the 
35-cm telescope at the University of New South Wales 
(UNSW) observatory in Sydney. These 3 polarimeters 
all use similar techniques, making use of ferro-electric 
liquid-crystal modulators operating at 500 Hz, and com- 
pact photomultiplier tubes (PMTs) as detectors. Some 
of the observations were previously reported by Cotton 
et al. (2016a) and Bailey et al. (2017), but all the data 
presented here have been re-analysed using the meth- 
ods described by Bailey et al. (2020b). We adopt a new 
calibration of modulator performance that has been car- 
ried out for an upcoming paper (Cotton et al. in prep.). 
'The resulting final polarization values are similar to past 
work. 

Details of the observing runs, made between 2014 and 
2023, as well as instrumental details are listed in Table 1. 
The polarization observations of e Sgr are listed in Ta- 
ble 2. New observations of interstellar control stars are 
described in Appendix A. The AAT observations were 
obtained with a range of broad-band filters. The Sloan 
Digital Sky Survey (SDSS) g’ and r' filters used with 
HIPPI were made by Omega Optics. For HIPPI-2 the 
corresponding filters were generation-2 filters from As- 
trodon Photometrics. The two Mini-HIPPI observations 
were made with no filter (Clear). A bandpass model as 
described in Bailey et al. (2020b) was used to determine 


the effective wavelength and modulation efficiency for 
each observation. This model takes account of all the 
optical components as well as the source spectrum, and 
atmospheric transmission. These values are listed in Ta- 
ble 2 and the effective wavelengths for these observations 
range from ~400 nm to ~720 nm. 


2.2. Imaging 


Infrared images of e Sgr were obtained on 2016 Sep 
10 using the IRIS2 instrument (Tinney et al. 2004) on 
the 3.9m Anglo-Australian Telescope in a narrow band 
filter at a wavelength of 2.3 um. The separation (p) 
and position angle of the companion were measured 
from the images with the Python Photutils package 
(Bradley et al. 2022) using profile fitting methods based 
on the DAOPHOT algorithms (Stetson 1987). The re- 
sults, together with those from past observations, are 
listed in Table 3. The separation of the companion has 
been slowly increasing, with little change to the position 
angle. 

'The companion was included in the aperture for all 
our polarization observations. However, it contributes 
less than 196 of the optical-region light, and as a main- 
sequence G star is not expected to have significant in- 
trinsic polarization (Cotton et al. 2017b). 


2.3. Archival Spectroscopy 


Spectra of e Sgr have been obtained from the archive 
of the Canada France Hawaii Telescope (CFHT). These 
were taken with the Echelle Spectropolarimetric Device 
for Observation of Stars (ESPaDOnS) instrument (Do- 
nati 2003) on 2014 Sep 8. The spectra cover 370 to 1000 
nm at R ~ 70000. The individual integrations have been 
combined to form a single spectrum which was used for 
the analysis in Section 3.3 and 3.4. 

We also used three spectra covering the Ha re- 
gion taken with the ultraviolet and visual echelle spec- 
trograph (UVES, Dekker et al. 2000) on the Euro- 
pean Southern Observatory (ESO) Very Large Telescope 
(VLT)!. These cover 472.6 to 683.5 nm and have R = 
14450. The dates of observation were 2016 Apr 22, May 
17 and Jun 1. 


3. MODELLING 
3.1. Polarization of Rotating Stars 


We model the polarization of rotating stars using the 
methods described in detail by Bailey et al. (2020a) and 
also used by Cotton et al. (2017a), Lewis et al. (2022) 


1 Based on data obtained from the ESO Science Archive Facility 


with DOI: 10.18727/archive/50 


Table 1. Summary of Observing Runs and telescope-polarization (TP) Calibrations 


Telescope and Instrument Set-Up^ Observations" Calibration® 
Run Date Range? Instr. Tel. f/ Ap. Mod. Filt. Det.^m Aeg Eff. qTP WTP 
(UT) (^) (nm) (X) (ppm) (ppm) 
2014AUG 2014-09-01 HIPPI AAT 8 6.6 BNS-E1 g 1463.1 89.0 —39.9 + 0.9 —38.2 + 0.9 


2016FEB/JUN^ 2016-02-27 to 06-25 HIPPI AAT 8 6.6 BNS-E2 g 2 463.8 86.9 —20.5 21.7 4541.8 


1619.9 806 -82+1.7 1341.7 


m2016JUN 2016-06-10 M-HIPPI UNSW 1158.9 MT Clear 1 467.9 75.5 —64.5 + 43 —7.6 + 4.4 
m2016JUL 2016-08-31 M-HIPPI UNSW 1158.9 MT Clear 1 464.0 74.3 —74.5 + 3.1 —11.7 + 3.1 
2017JUN 2017-06-26 to 07-04 HIPPI AAT 8 6.6 BNS-E2 425SP 1400.1 49.4 -764+3.8 88+ 3.8 
9g 1 463.3 86.8 —9.1 + 1.5 —2.6 + 1.4 

650LP 1717.0 62.7 —8.6 + 2.5 —5.3 + 2.5 

2017AUG 2017-08-08 to 08-14 HIPPI AAT 8 6.6 BNS-E2 425SP 1400.2 49.5 -764+3.8 88:338 


500SP 1 435.3 74.0 —10.0 + 1.7 —0.4 + 1.6 
2619.2 80.7 —10.6 + 13 —7.1 + 1.3 
2 461.7 81.5 128.1 +0.9 3.8 + 0.9 
2620.9 84.5 109.1 + 1.4 6.91.4 
1532.3 96.6 —20.1 + 1.5 2.3 1.5 
1602.0 88.9 —10.3 + 2.2 3.6 + 2.2 


1 438.2 55.1 1.9 1.4 18.0 + 1.4 


2018MAR 2018-03-18 to 04-06 HIPPI-2 AAT 8* 25.5 BNS-E3 g 


2018JUL? 2018-07-12 to 07-25 HIPPI-2 AAT 8*11.9 BNS-E4 V 


2018AUG? 2018-08-21 to 09-01 HIPPI-2 AAT 8* 11.9 BNS-E5 500SP 


V 1532.4 95.7 —20.1 + 1.5 2341.5 

650LP 1720.9 76.8 -65+1.8 38-18 

2023APR/MAY^ 2023-06-03 HIPPI-2 AAT 1512.7 ML-EI 425SP 1398.8 79.5 -1042.4 7.8 2.1 
500SP 1 438.2 86.9 —20.1 + 1.7 —5.1 + 1.5 


1460.2 89.2 5.9+1.0 0.0 + 1.0 
1 601.8 64.2 4.9 + 4.5 —19.7 + 3.9 
1 620.8 60.5 —2.7 + 3.2 —18.8 + 3.3 
1719.1 44.6 —36.6 + 5.8 —35.4 + 6.2 


E» 
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650LP 


Notes: 

* Indicates use of a 2x negative achromatic lens, effectively making the foci f/16. 

^ A full description, along with transmission curves for all the components and modulation characterisation of each modulator 
(‘Mod.’) in the specified performance era, can be found in Bailey et al. (2020b). 

^ Mean values are given as representative of the observations made of e Sgr. Individual values are given in Table 2 for each 
observation; n is the number of observations of e Sgr. 

^ Most of the observations used to determine the TP and the high-polarization standards observed to calibrate position angle 
(PA), are described in Bailey et al. (2017) (m2016JUN, m2016JUL), Bailey et al. (2015) (2014M AY), Cotton et al. (2017b) 
(2016FEB, 2016JUN) Cotton et al. (2019) (2017JUN, 2017AUG), and Bailey et al. (2020b) (other runs). This is the first data 
presented from the 2022M AY run; 2-6 total observations per band of HD 2151, HD 61421, HD 102647, HD 102870 or HD 
140573 were used to determine TP; PA was calibrated with reference to HD 84810, HD 161471, HD 111613 and HD 183143. 

1 Dates given are for observations of e Sgr and/or control stars. 

* B, R indicate blue- and red-sensitive H10720-210 and H10720-20 photomultiplier-tube detectors, respectively. 

f TP calibration was carried out using observations combined from 2016FEB and 2016JUN runs. 

? TP calibration was carried out using observations combined from 2018JUL and 2018AUG runs. 

^ TP calibration was carried out using observations combined from 2023APR A/B and 2023MAY runs. 


Table 2. Polarization Observations of e Sgr (Sorted by Effective Wavelength). 


Run UT (mid-point) Dwell Exp. Filt. Det.^ Acq Eff. q u p 0 
(s) (s) (nm) (%) — (ppm) (ppm) (ppm) (°) 
2023MAY 2023-06-03 16:39:20 1615 1280 425SP B 398.8 79.3 —5.8 + 12.1 —13.6 + 12.0 14.8 + 12.0 123.6 + 27.9 
2017JUN 2017-06-26 16:29:55 1557 800 425SP B 400.1 49.4 6.7 + 15.6 —50.3 + 15.6 50.7 + 15.6 138.8 + 9.2 
2017AUG 2017-08-10 14:14:07 2124 800 425SP B 400.2 49.5 28.1 + 14.6 —19.4 + 14.5 34.1 + 14.6 162.7 + 14.2 
2023MAY 2023-06-03 17:03:20 996 640 500SP B 433.4 86.9 36.0+ 6.7 600.40 6.7 70.3+ 6.7 29.6 + 2.7 
2017AUG 2017-08-11 13:53:51 1506 640 500SP B 435.3 74.0 21.8 + 8.9 810+ 89 839+ 89 37.5+ 3.0 
2018AUG 2018-08-21 14:16:31 973 640 500SP B 438.2 54.9 25.1+ 6.8 83.2+ 6.8 86.9+ 6.8 36.6+ 2.2 
2023MAY 2023-06-03 16:17:39 983 640 g B 460.2 89.2 19.304 2.5 95.1+ 2.6 9704 2.6 39.3+ 0.8 
2018MAR? 2018-04-06 17:39:10 1032 640 g B 461.7 81.5 71.34 2.5 99.9404 2.5 122.7 2.6 27.24 0.6 
2018MAR 2018-04-06 17:19:21 1189 640 g' B 461.8 81.6 162+ 2.9 955+ 2.7 969+ 2.8 40.2 0.8 
2014AUG 2014-09-01 10:11:12 1843 640 g B 463.1 89.0 41.04 3.6 164.0 + 3.6 169.0 + 3.6 38.0 0.6 
2017JUN 2017-06-26 16:02:50 1588 800 g B 463.3 86.8 13.4+ 6.0 102.0 + 4.2 1029+ 5.1 41.3 1.7 
2016FEB 2016-03-01 18:37:33 973 480 g B 463.8 86.9 276+ 3.7 888+ 3.7 93.0 + 3.7 36.4 1.2 
2016FEB 2016-02-29 18:34:56 1012 480 g B 463.9 87.0 36.2+ 4.0 851+ 4.1 9254 4.0 335+ 1.2 
m2016JUL 2016-08-31 10:37:32 1609 800 Clear B 464.0 74.3 43.5 19.9 123.2 + 20.6 130.7 + 20.3 35.34 4.4 
m2016JUN 2016-06-10 10:47:32 1850 800 Clear B 467.9 75.5 29.5 + 21.9 98.45 + 21.8 102.7 + 21.8 36.74 6.2 
2018JUL 2018-07-25 16:08:13 1025 640 V B 532.3 96.6 39.8 + 4.1 217.1 + 4.1 220.074 4.1 39.8+ 0.5 
2018AUG 2018-08-21 14:36:42 1345 960 V B 532.4 95.7 45.204 3.1 178.5 + 3.0 184.1 + 3.1 37.9 0.5 
2023MAY 2023-06-03 17:23:21 1405 960 r B 601.8 64.2 31.5 + 7.5 261.8 + 6.7 263.74 7.1 41.6 0.8 
2018JUL 2018-07-25 16:26:18 1045 710 r B 602.0 88.9 51.04 7.3 262.8 + 6.4 267.7 + 6.8 39.5 0.8 
2017AUG 2017-08-09 13:12:24 1880 640 r R 619.2 80.7 61.1 + 3.9 281.3 + 3.7 287.9+ 3.8 38.9+ 0.4 
2017AUG 2017-08-08 14:18:24 1339 640 r R 619.3 80.7 63.8 + 3.7 283.9 + 3.7 291.0 + 3.7 38.7 t 0.4 
2016FEB 2016-02-27 16:41:49 1237 640 r R 619.9 80.6 30.3 + 5.2 227.5 + 5.4 229.5 + 5.3 41.20 0.7 
2023MAY 2023-06-03 18:24:03 958 640 r R 620.8 60.5 645 + 5.1 284.8 + 5.2 292.0 + 5.2 386+ 0.5 
2018MAR! 2018-03-18 18:49:50 980 640 r R 620.9 84.5 784-4 3.2 276.1 + 3.2 287.0 + 3.2 37.1+ 0.3 
2018MAR? 2018-03-18 19:11:58 1300 640 r R 620.9 84.5 71.94 4.6 278.7 + 4.4 287.8 + 4.4 37.8 0.5 
2017JUN 2017-07-04 13:08:09 2165 1440 650LP R 717.6 62.7 81.9 4.8 3285+ 4.8 3386+ 4.8 380+ 0.4 
2023MAY 2023-06-03 18:45:19 1634 1280 650LP R 719.1 44.6 101.34 8.1 341.24 8.5 35594 8.3 36.74 0.7 
2018AUG 2018-08-23 11:47:56 1232 960 650LP R 720.9 76.8 524+ 4.1 318.24 4.1 322.5 + 4.1 40.34 0.4 


^ B, R indicate blue- and red-sensitive H10720-210 and H10720-20 photomultiplier-tube detectors, respectively. 


s/l Aperture 


sizes are as given in Table 1 except where indicated by these superscripts; l: 25.5 arc seconds, s: 5.3 arc seconds — small 
enough to exclude the reported debris disk but not the binary companion (Rodriguez & Zuckerman 2012). 


Table 3. Measurements of the e Sgr Binary System. 


UT Date — p(") pa (°) Am (mag) Ref 

1992-06-18 2.06+0.01 142+1 5.85 (V) 1 

1999-03^ — 2.39 142 — 4.6 (K) 2 

2016-09-10 2.58+0.04 140+1 4.2 (2.3um) 3 
Notes: 


References: 1 Golimowski et al. (1993), 2 Hubrig et al. 
(2001), 3. This paper (Section 2.2) 
* The exact date of the observation is not reported. 


and Howarth et al. (2023). We use a Roche model for the 
rotating star and assume gravity darkening according 


to Espinosa Lara & Rieutord (2011). This allows us 
to determine the distribution of local temperature and 
gravity as a function of latitude. We use a set of custom 
ATLASO (Castelli & Kurucz 2003) models to represent 
the local stellar atmosphere at 46 colatitudes from 0 to 
90? at 2? steps. For each of these model atmospheres we 
then calculate the specific intensity and polarization of 
the emergent radiation. These calculations use a version 
of the SYNSPEC spectral synthesis code (Hubeny et al. 
1985; Hubeny 2012) which we have modified to do full 
polarized radiative transfer using the vector linearized 
discrete ordinate radiative transfer (VLIDORT) code of 
Spurr (2006). 


Eps Sgr polarization (w = 0.997, inc = 90°) 
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Figure 1. Polarization observations (rotated to put all the 
polarization into the Q stokes parameter) compared with a 
model of polarization for a rapidly rotating star. The ob- 
servations show similar wavelength dependence to the model 
but the polarization is offset by ~300 ppm in the positive 
direction. 


We overlay a rectangular grid of pixels over the pro- 
jected view of the star. For each pixel we determine the 
local temperature, gravity, and the surface normal view- 
ing angle. We then use these values to interpolate in our 
set of radiative transfer calculations to determine the lo- 
cal specific intensity and polarization. Summing over all 
the pixels then gives the integrated flux and polarization 
as a function of wavelength. 


3.2. Polarization Components 


Figure 1 shows the polarization observations rotated 
through ~40°, which puts the polarization primarily in 
the Q Stokes parameter. The polarization shows a large 
change with wavelength from values of about -50 ppm at 
7-400 nm, to more than 300 ppm at ~720 nm. This large 
wavelength change matches that expected for a star with 
a large rotation rate. However, when we compare the ob- 
servations with a rapidly-rotating-star model, as shown 
in Figure 1, we find that the observed polarization is 
offset in a positive direction from the modelled value by 
~300 ppm. In this plot negative polarization is perpen- 
dicular to the star's rotation axis, positive polarization 
is parallel to the axis. 

For the polarization to be consistent with rotational 
distortion, there must be an additional source of polar- 
ization of ~300 ppm with a position angle parallel to 
the star’s rotation axis. In our previous studies of rotat- 
ing stars we have been able to model the observations 
using a rotating star model and a small contribution of 
interstellar polarization. For e Sgr at a distance of 44 pc 
(van Leeuwen 2007) the expected level of interstellar po- 
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larization is ~60 ppm according to the models of Cotton 
et al. (2017b) and a value as high as 300 ppm would be 
extremely unlikely (see Appendix A). Furthermore the 
interstellar polarization would have to be aligned with 
the rotation axis of the star to provide the required off- 
set. This would not be consistent with what is measured 
for nearby stars and seems intrinsically unlikely. 

For e Sgr there is, however, an additional likely source 
of polarization with the required properties. The narrow 
shell absorption features in the spectrum, as described in 
Section 1, indicates the presence of circumstellar gas. In 
such a rapidly rotating star this is most likely in the form 
of a narrow equatorial disk. Such disks are common in 
rapidly rotating B stars as seen in the classical Be stars 
(Rivinius et al. 2013). One of the important observable 
features of Be stars is polarization due to starlight scat- 
tered from the disk. The typical polarization level in 
nearby Be stars (where interstellar polarization is unim- 
portant) is ~3000-6000 ppm (Cotton et al. 2016a). 

The much lower polarization (7-300 ppm) we infer for 
the disk in € Sgr, as well as the lack of strong Ha emis- 
sion, indicates a disk of much lower density than those 
usually found in classical Be stars. At low densities we 
expect electron scattering to be the dominant opacity 
source leading to a flat, wavelength-independent polar- 
ization (Halonen & Jones 2013; Rivinius et al. 2013). 
This differs from the wavelength-dependent polarization 
structure seen in the higher-density disks of classical 
Be stars (Poeckert, Bastien, & Landstreet 1979; Wood, 
Bjorkman, & Bjorkman 1997) where bound-free absorp- 
tion is important. Comparison with the models of Halo- 
nen & Jones (2013) suggests a gas density of ~ 1071? g 
cm? whereas densities up to ~ 10710 g cm^? are seen 
in classical Be stars. 

'The polarization position angle of a Be star disk is per- 
pendicular to the disk (Quirrenbach et al. 1997). Hence 
the polarization is parallel to the star's rotation axis, as 
required for our offsetting polarization component. 


3.3. Polarization Variability 


The polarization observations of e Sgr (Table 2) show 
considerable scatter between observations in the same 
filter on different dates. The total range in p values is 
77 ppm at g' and 62 ppm at r’. This is much larger than 
the statistical errors, and much larger than the scatter 
seen in our similar observations of other rapidly rotat- 
ing stars using the same instruments (e.g. Dailey et al. 
2020a; Howarth et al. 2023). The most extreme values 
are seen in some of the earliest observations in 2014 and 
2016, whereas observations closely spaced in time gen- 
erally show good agreement. This suggests variability 
of the polarization on long timescales. Such variability 
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Figure 2. Spectra of e Sgr in the Ha region. The UVES 
(ESO VLT) spectrum is from 2016 Apr 22 and shows broad 
absorption and the narrow shell absorption component. The 
ESPaDOnS (CFHT) spectrum on 2014 Sep 8 shows addi- 
tionally double-peaked core emission. The very narrow ab- 
sorption features are telluric lines. 


is not unexpected for a system with a disk component 
as described in the previous section. Polarization in Be 
stars is often variable (Coyne 1976). 

Independent evidence for a variable Be-type disk in € 
Sgr comes from the archival spectroscopy (Section 2.3). 
Figure 2 shows two spectra of e Sgr. One of these, the 
CFHT spectrum from 2014 Sep 8, shows weak double- 
peaked emission in Ha. The other spectrum, taken 
on 2016 Apr 22 with VLT UVES shows a broad pho- 
tospheric absorption and narrow shell-absorption core. 
Similar behavior is seen in the other UVES spectra from 
2016. Since the earlier spectroscopy discussed in Sec- 
tion 1 makes no mention of emission we assume the 2016 
spectra represent the more usual state. Interestingly the 
highest g^ polarization was observed on 2014 Sep 1, only 
a few days before the spectrum showing the emission, 
whereas polarization observations in 2016 show some of 
the lowest values. 


3.4. Model Grid 


The polarization of a rotating-star model (as described 
in Section 3.1) depends on four main parameters: the 
polar temperature and gravity (Tp, gp), the rotation rate 
(w = Q/Qcrit)?, and the inclination of the star's rotation 
axis (i). The polarization observations alone are not 
sufficient to determine these four parameters so, as in 


2 Where the subscript ‘crit’ is used to denote critical rotation 


throughout this paper. 
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Figure 3. Results of Fourier- Transform modelling of CFHT 
spectra of the MgII 448.1nm line to determine ve sini. The 
value plotted is the rms difference between observed and 
modelled Fourier Transforms. 


our previous studies (Bailey et al. 2020a; Lewis et al. 
2022; Howarth et al. 2023) we use additional information 
to provide relationships between these parameters. 

'The first of these is the spectral-energy distribution, 
which principally constrains the temperature and an- 
gular diameter (i.e., radius). We use photometry from 
Johnson et al. (1966) (optical) and archival spectroscopy 
obtained with the International Ultraviolet Explorer 
(IUE; UV). IUE spectra taken in low-resolution mode 
(AA ~ 0.6 nm) through the spectrographs’ large aper- 
ture provide the most reliable spectrophotometry. Only 
one spectrum is available for each of the short- and long- 
wavelength ranges (AA ~ 115-198, 185-335 nm), but the 
data are well exposed (SWP 16443, LWR 12689). We 
adjusted the fluxes to match the current HST calibration 
by using results from Massa & Fitzpatrick (2000), in- 
cluding the minor update described by Fitzpatrick et al. 
(2019). This gives an ~8% increase in integrated UV 
flux over the standard pipeline product (and a ~1.5% 
increase in the inferred Tyg). The resulting integrated 
observed flux in the adopted 123.5-321 nm region is 1.40 
x 10-9 erg cm7? s71. We then applied an extinction 
correction for E(B — V) = 0008, estimated using the 
G-TOMO tool? (based on results from Lallement et al. 
2022), with a Seaton (1979) extinction curve (adopting 
Rz A(V)/E(B — V) = 3.1). 

The second constraint is the value of vesini. We de- 
termined new values of vesini using a Fourier Trans- 
form analysis of archival CFHT spectra. We used the 
MglI 448.1 nm line, chosen as the rotational line shape 
is less affected by shell absorption than alternatives such 
as Sill 634.7 nm. The analysis was done by comparing 
model spectra, calculated using Roche geometry, to the 


3 https://explore-platform.eu/sdas 
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Figure 4. Fit to e Sgr grid using the rotational polarization 
models plus a fixed polarization offset. In this and similar 
figures each point on the grid corresponds to a set of stel- 
lar parameters that fits the spectral energy distribution and 
ve sin i of e Sgr as described in Section 3.4. The plot shows 
how well each of these grid models fits the polarization ob- 
servations. 


observations in Fourier-Transform space (see Howarth 
et al. 2023, for a more detailed description of this ap- 
proach). Because of the rotational distortion of the star 
and the strong gravity darkening the inferred ve sin? val- 
ues depend on the rotation rate w, as shown in Figure 3. 
For slow rotations (alone) the inferred v, sin i is in rea- 
sonable agreement with the value of 236 km s^! given 
by Royer et al. (2002). 

We adopt the ve sin iw relation of Figure 3, together 
with the the Hipparcos parallax of 22.76-50.24 mas 
(van Leeuwen 2007; the star is too bright for Gaia). 
We can then construct a two-dimensional grid of stel- 
lar parameters where all other values can be established 
from the two parameters w and i. Each model in this 
grid corresponds to a star that fits the spectral-energy 
distribution and vesini. We can then calculate the po- 
larization for each of these grid models and compare 
with the observations. 


3.5. Comparison with Observations 


In comparing the model grid calculations with the ob- 
servations we have to contend with the variability of 
the system (Section 3.3) and the presence of additional 
polarization components due to the disk and interstel- 
lar polarization (Section 3.2). There is insufficient data 
to attempt to correct for the variability. However, ex- 
amination of the polarization measurements shows that 
the earliest observations (those in 2014 to 2016) seem to 
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Figure 5. As Figure 4 but using the rotational polarization 
models plus a disk and up to 30ppm interstellar component 
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Figure 6. As Figure 4 but using the rotational polarization 
models plus a disk and up to 100 ppm interstellar compo- 
nent. The solid lines show the 2c contours on the rotational 
parameters as fitted to these models (see Section 3.6). The 
dash-dot lines shows the corresponding extended range of pa- 
rameters when the uncertainties on the input parameters to 
the model grid (in particular ve sin i) are taken into account. 


show the largest scatter, and variability at these times is 
confirmed by the spectra shown in Figure 2. We there- 
fore chose to exclude these observations from the analy- 
sis, which was restricted to the 22 observations covering 
2017 — 2023. There is still some evidence of variabil- 
ity, which shows up as excess noise compared with that 
expected from the measurement errors. 
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Table 4. Parameters of Best Fit Region in Grid Model 
Variants (Figures 4 to 6) 


Variant QU Offset Disk + Interstellar 

Pmax (ppm) 0.0 30.0 60.0 100.0 

lism (deg) 144.3 137.1 134.2 

pais (ppm)^ 276 299 330 367 

brot (deg) 40.1 39.9 39.9 39.9 

w(i = 90) 0.9977(7) 0.9975(7) 0.9972(7) 0.9968(9) 

w(i = 88) 0.9977(7) 0.9974(7) 0.9971(7) 0.9967(9) 

w(i = 86) 0.9981(6) 0.9978(7) 0.9975(7) 0.9971(8) 

w(i = 84) 0.9986(6) 0.9983(7) 0.9979(6) 0.9974(7) 

w(i = 82) 0.999 0.999 0.9987(5) 0.9982(6) 

w(i = 80) 0.999 0.999 
Notes: 


Figures in parentheses are the errors in the last digit of the 
value. When the best fit for an inclination is in the 
w=0.999 bin the value give is the bin centre and no error 
can be determined. 

^ Varies slightly with inclination, average value is given. 


Table 5. Best-fitting Stellar Parameters. 


Parameter i=90° Best? Range’ (20) 
w 0.997 0.999 0.995 — 1.0 
Ve/Verit 0.95 0.97 0.94 — 1.0 
Inclination, i [°] 90 80 76 — 90 
Ter [K] 10091 9950 9845 — 10096 
Tp [K] 11791 11721 11548 — 11813 
T. [K] 7884 7433 7324 — 1492 
Log(L/Lo) 2.709 2.696 2.669 — 2.726 
R, [Rol 5.98 6.01 5.90 — 6.08 
Re [Ro] 8.59 8.80 8.47 — 8.89 
Mass [Mo] 3.69 — 3.80 3.69 — 3.90 
log g» [dex cgs] 3.45 3.46 3.45 — 3.48 
log ge [dex cgs] 2.24 2.00 1.99 — 2.36 
P; [days] 163 161 1.57 — 1.64 
Notes: 


“ The parameters listed are those of the grid model (see 

Section 3.4) giving the lowest x? among all models (Best), 

and the best model at inclination of 90 degrees. 

oR ; : 
ange is the full range of stellar parameters corresponding 

to the region enclosed by the dot-dash line in Figure 6. 


'The additional polarization components were handled 
by fitting the observations with several variants of the 
model. The same grid of calculations was used for the 
polarization of the distorted star, but different models 
were used for the polarization of the additional com- 


ponents. Initial tests over a wider and coarser grid of 
models showed that only the models with high w were 
compatible with the polarization. The final grid covered 
the range 0.980 < w < 0.999 and 70 < i < 90 and the 
relevant section of this region is shown in Figures 4 to 
6. 

In the first variant (Figure 4) we represent the addi- 
tional polarization as a fixed (i.e. wavelength indepen- 
dent) offset in Q and U. Thus at each point in the grid 
we fit three values, the Q-offset, U-offset and the posi- 
tion angle of the star’s rotation axis (rot). With this 
model we find that for the best-fitting (lowest-X?) points 
in Figure 4, the additional QU-offset corresponds to a 
position angle of 38.5°, and the star’s rotation axis is at 
40.1°. The polarization is within 1.6° of that expected 
for an equatorial disk (Section 3.2). The small difference 
in angles can be understood as due to an additional con- 
tribution from interstellar polarization. 

In the remaining variants we represent the addi- 
tional polarization as a combination of a wavelength- 
independent disk contribution, and an interstellar- 
polarization component following the ‘Serkowski law’ 
(Serkowski 1973; Serkowski et al. 1975; Wilking et al. 
1982; Whittet et al. 1992). We set the Amax of the inter- 
stellar polarization to 470 nm, the typical value found 
within the Local Hot Bubble (Marshall et al. 2020 and 
refs. within). In these models there are four fitted pa- 
rameters at each grid point: pmax and the position angle 
Jism of the interstellar polarization, the position angle 
Orot of the star’s rotation axis, and the polarization of 
the disk paisk. The disk is assumed to be in the equato- 
rial plane so its polarization position angle is 0... 

There is some degeneracy in these models between 
Pmax and paix. We therefore ran three versions of 
the models where we applied the restricted ranges 0 « 
Pmax € 30 ppm, 0 < pmax < 60 ppm, and 0 < pmax < 
100 ppm, based on our analysis of the likely value of in- 
terstellar polarization as described in Appendix A. Al- 
though we constrain pmax to a range, the best-fitting 
models in all these cases have Pmax at the largest allowed 
value. The results for two of these cases are plotted in 
Figures 5 and 6, and all the variants are described in 
Table 4. 

All these variants (Figures 4-6) produce quite similar 
results, with a “valley” of low x? values running from 
w ~ 0.997 at i = 90° to w = 0.999 at à ~ 80°. Thus 
despite different assumptions about the interstellar po- 
larization, almost the same stellar parameters are being 
selected in all these cases. 

In Table 5 we list the stellar parameters corresponding 
to the two grid models at the ends of this x? valley. The 
model labelled *Best" (in the *Disk 4- 100 ppm' set) is 
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Figure 7. Polarization and intensity images for the two 
models of e Sgr listed in Table 5. Specific intensity is overlaid 
with polarization vectors for a wavelength of 400 nm. 


the one that has the lowest x? and is at w = 0.999, i = 
80° (in some of the other model variants the best fit is in 
the adjacent grid point at i = 82°). However, the very 
strong shell absorption seen in e Sgr, despite the low 
density of the disk, might favour a higher-inclination 
model like the w = 0.997, i = 90° case also listed in 
Table 5. 

Figure 7 shows images of the stellar models for the two 
cases listed in Table 5. The models show the distribution 
of specific intensity over the star, showing the strong 
effect of gravity darkening at the equator, and is overlaid 
with polarization vectors. The polarization is dominated 
by light from the bright regions at top and bottom of 
the image resulting in the negative (perpendicular to the 
rotation axis) integrated polarization at this wavelength. 


3.6. Uncertainties 


To determine the uncertainties on the fitted values 
of w and 4 we use the bootstrap method (Press et al. 
2007). For each of the model grid variants we ran the 
modelling on 1000 trials using random selections of the 
22 observations with replacement (so that each obser- 
vations can be selected more than once). For each trial 
we determined the best w value at each inclination by 
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fitting the x? values with a cubic function. We then de- 
termined the standard deviation of these values over the 
1000 trials. These results give the errors on w listed in 
Table 4 and plotted as error bars on Figure 4 to 6. This 
method of error determination does not depend on the 
errors on the individual polarization measurements, but 
on the actual scatter on the data, which is largely due 
to the polarization variability of the star. 

Using these bootstrap results we can determine the 
Ax? ranges corresponding to 1 or 2c, and thus deter- 
mine the 2-c contours on w and 7 which are shown as 
the solid lines in Figure 6. We use Figure 6 (the ‘Disk + 
100 ppm’ variant) as this gives the lowest x? values, and 
has a somewhat broader distribution in w, but largely 
overlays the best-fit regions covered by the other model 
variants. 

Additional uncertainties are introduced by the uncer- 
tainties in the input parameters to the model grid (Sec- 
tion 3.4). These were assessed by running the mod- 
elling for auxiliary grids calculated for adjusted values 
of ve sini, parallax and UV flux in turn. The most sig- 
nificant effect is that due to an estimated 10 km s^! 
uncertainty in ve sin ? which has the effect of moving the 
x? valley up and down in w. This leads to the wider 2-o 
contours shown by the dash-dot lines in Figure 6. Un- 
certainties in parallax and UV flux do not significantly 
change the fitted values of w and i, but do have a small 
effect on the relationship to other stellar parameters. 
The wider contours (dash-dot line) were then used to 
determine the allowed range of each stellar parameter, 
as listed in Table 5. 


4. DISCUSSION 
4.]. There is No Polarizing Debris Disk 


Previously, HIPPI observations of e Sgr were reported 
by Cotton et al. (2016b). They interpreted the observed 
polarization in the context of the system's reported in- 
frared excess (Rodriguez & Zuckerman 2012). Cotton 
et al. present a model comprising two dust components 
(to reproduce the continuum IR excess) in addition to 
the stars. A cold dust belt (responsible for the bulk of 
the IR excess) was assumed to be circumbinary, with a 
warm dust belt orbiting one of the two stars. A residual 
polarization of 30 ppm remained after accounting for po- 
larization induced by the interstellar medium. To repro- 
duce this a fractional polarization of around 5096 of the 
total IR excess would be required. Such a value is much 
higher than those measured for dust grains in other de- 
bris disks, which peak at around 20 to 30% (e.g. Krivova 
et al. 2000; Graham et al. 2007; Anche et al. 2023). Since 
the Cotton et al. (2016b) model was proposed, work by 
Vandeportal et al. (2019) has shown that warm (mid-IR) 
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excesses represent substantially lower dust masses and 2300 ppm, compared with «100 ppm for the other 
contribute less to polarization compared to cool (far- rapidly-rotating stars (Cotton et al. 20172; Bailey et al. 
IR) excesses with lower fractional luminosities but rep- 2020a; Howarth et al. 2023). The only other star with 
resenting greater dust masses. This further reduces the a comparable polarization amplitude was 0 Sco (Lewis 
plausibility of the proposed disk model. et al. 2022). In that case the rotation was slower than in 
The large beam size of IRAS and the low galactic lat- € Sgr, but the low gravity of this evolved star enhanced 
itude of the system (|b| < 10°) also lead to the possi- the polarization. 
bility that the flux density measurement at 60 um is Multiwavelength polarimetry should therefore be a 
contaminated by diffuse background emission. Indeed, good way of finding other similar extreme rotators, and 
the IRAS point source catalog (Joint Iras Science 1994) establishing how common such objects are. In the 19 
notes a quality flag of “2” for the 60 um band. Applying stars we have observed for our AAT polarimetry pro- 
a downward revision of the flux, as implied, would only gram on rotating stars we have found no other candi- 
increase the tension between the level of IR excess and dates for such extreme rotation. The program is based 
the polarization efficiency of the disk dust required. on stars down to V ~ 3.5 from the list of van Belle 
In the light of observations presented here we can (2012). 
eliminate the disk model as a plausible explanation for 
the polarimetric signal arising from the system. Specif- 4.3. Relationship to Be Stars 


ically, the invariance of the position angle of the polari- 
metric signal between epochs, despite the secondary's 
proper motion, is inconsistent with illumination of a cir- 
cumbinary dust disk. The agreement between the r’ 
observations with small and large aperture sizes in the 
2018MAR run is also inconsistent with this hypothesis. 
Furthermore, the level of fractional polarization required 
would be a factor of two higher than other systems for 
which such measurements are available (see above, and 
also Marshall et al. 2023). 


Given that Rodriguez & Zuckerman (2012) find the than this, with a broad distribution of ve/ Verit values 
from ~0.5 to ~0.95, with the mode of the distribution at 


0.66 (Zorec et al. 2016). There are also rapidly-rotating 
B stars (e.g. Regulus, Cotton et al. 2017a) that are not 
Be stars. 

While e Sgr has a rotation rate at the top end of the 
range inferred for Be stars it is certainly not a typical 
classical Be star. It does appear to have an equato- 
rial disk, but it is very weak compared with typical Be 
4.2. Rapid Rotation stars. Weak Ha emission has been observed on one oc- 
casion but is normally absent, as discussed in section 
3.3. Thus, these results confirm that rotation alone can- 
not be the key requirement for formation of Be disks. 
Rather, there must be some additional mechanism op- 
erating in Be stars that can launch material from the 
surface into orbit. In e Sgr this mechanism is either not 
operating or is operating with much reduced efficiency, 
enabling the formation of only a weak disk despite the 
very rapid rotation. We note that Classical Be stars are 
more common at early B spectral types and less common 
at the B9/A0 spectral type of e Sgr. Huang et al. (2010) 
suggest that the threshold rotation rate for forming a Be 
disk is much higher for later type stars. 

Most classical Be stars show photometric variability 
^ VFTS285, an O star in the Large Magellanic Cloud, has been on short timescales. Observations with the Transiting 

reported as having w = 1.00 (Shepard et al. 2022) Exoplanet Survey Satellite (TESS, Ricker et al. 2015) 


The value of ve/Verit ~0.95 that we find for e Sgr is 
interesting because this is the value that was suggested 
(Townsend et al. 2004; Owocki 2005) as required for 
“weak” processes such as pulsation to be able to launch 
material from the star into orbit. If all Be stars rotated 
this fast, the existence of Be disks would be easier to ex- 
plain. However, statistical analysis of v, sin i data for Be 
stars (Cranmer 2005 [but see Howarth 2007]; Zorec et al. 
2016) suggests that most Be stars rotate much slower 


system to be unstable with the accepted debris disk 
architecture and that IR. excess may also be produced 
by free-free emission in the (previously unidentified) Be 
disk (Kastner & Mazzali 1989), it is likely no debris disk 
exists at all. The nature of e Sgr as a rapidly rotating 
B star with a gas disk is therefore the origin of its po- 
larization. 


'The rotation as a fraction of criticality that we mea- 
sure for e Sgr of w ~ 0.997, which corresponds to ve / vcri 
70.95, is the largest yet measured for a star in our 
galaxy*. The previous highest was obtained by Domi- 
ciano de Souza et al. (2014) who measured the rotation 
of Achernar using interferometry as w = 0.980, Ve/Verit 
= 0.883. Achernar is a Be star, but the measurements 
were made during a quiescent phase where the Be disk 
was absent. 

€ Sgr stood out in our polarization study of bright, 
rapidly-rotating stars because of the large amplitude 
of the wavelength-dependent polarization variations — 


11 


€ Sgr polarization (inclination = 90°) 
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Figure 8. Comparison of observations and models for the polarization of e Sgr as a function of wavelength. The models 
show the predictions due to stellar rotation. The disk and interstellar polarization have been applied as a correction to the 
observations. The black solid line is the model for w — 0.997, i — 90?. Blue diamonds are the predicted polarization values for 
each observed filter (they are not necessarily exactly on the line, because they are averages over the filter bandpass). Grey lines 


show the model predictions for grid models with lower values of w. 


typically show multiple periodic signals at frequencies 
between 0.5 and 4 day"! (Labadie-Bartz et al. 2022). 
These are normally interpreted as non-radial pulsations, 
which may be related to the mass ejection process that 
forms the disk (Rivinius et al. 2013). 

In Figure 9 we show the periodogram of the TESS two- 
minute cadence light curve of e Sgr for sector 13 (2019 
Jun 19 — Jul 17)°. The pattern seen of two groups of 
peaks (here at ~0.65 and 1.3 days~'), with a frequency 
ratio of roughly 1 to 2 is one that is very commonly seen 
in Be stars (e.g. Walker et al. 2005; Semaan et al. 2018; 
Baade et al. 2018; Labadie-Bartz et al. 2022). However 
the amplitude (around 140 ppm for the highest peak) is 
low compared with those more typically seen in classical 
Be stars, which are in the parts per thousands (Labadie- 
Bartz et al. 2022). 


5 The TESS data used here are available at MAST: 10.17909/t1ae- 


rd45 


As an example we show in the lower panel for Figure 
9 the periodogram of e Sgr compared with that of the 
bright Be star Achernar (aœ Eri, HD 10144), based on 
observations obtained from TESS sector 2 (2018 Aug 
23 — Sep 20). At that time Achernar was in an active 
state. It has been quiescent without a disk since late 
20199. The frequency structure seen in Achernar is quite 
similar to that of e Sgr, but the peak amplitude is 25 
times higher in Achernar. 


4.4. Evolutionary Status 


In Figure 10 we show the position of e Sgr on an HR 
diagram compared with evolutionary models for rotating 
stars of solar metallicity (Z — 0.014) from Georgy et al. 
(2013). The models are for M/Mg of 3.0 and 4.0 (our 
measured value is M/Mọ = 3.69), and initial rotation 


6 According to our polarization monitoring of Achernar over the 


last 10 years. 
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Figure 9. Periodogram of TESS light curve of e Sgr (TESS 
sector 13). In the upper panel it is shown at full scale. In 
the lower panel it is compared with the periodogram of the 
bright classical Be star Achernar (a Eri) from TESS sector 2. 
'The presence of two frequency groups with a 2-to-1 frequency 
ratio is common in Be stars, but the amplitude seen in € Sgr 
is much lower than is commonly seen. 
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Figure 10. Location of e Sgr on the HR diagram compared 
with evolutionary models for rotating stars from Georgy 
et al. (2013). Isochrones are shown at 10 Myr intervals from 
150 to 190 Myr, and then at 196 and 197 Myr. 
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Figure 11. Plot of w against log(g,) for e Sgr and the same 
models as in Figure 10. 'The rotation rate (w) reaches a 
maximum at a point corresponding to the “hook” in the HR 
diagram, and then declines rapidly as the star evolves across 
the Hertzprung gap. 


rates of w = 0.7, 0.8 and 0.9. The location of e Sgr is that 
of a star just beginning to evolve away from the main 
sequence. Its position corresponds to the hydrogen shell 
burning phase where the star is crossing the Hertzprung 
gap. This is a very rapid phase of stellar evolution and 
such stars are not commonly observed. 

Figure 11 shows the comparison with the same set of 
models in a plot of w against log(g,). The log(g,) we 
measure for e Sgr is again consistent with a star begin- 
ning to evolve away from the main sequence. The max- 
imum rotation rate for these model stars corresponds to 
the *hook" in the HR diagram where the star contracts 
after the end of core hydrogen burning. At this point 
a star born with w — 0.9 or greater reaches near criti- 
cal rotation. However, Figures 10 and 11 show that our 
measured parameters for e Sgr place the star at a point 
beyond the “hook” where the star is expanding and the 
rotation rate is slowing. The rotation we measure for 
€ Sgr is therefore too high to match any of these single 
rotating star models. Its rapid rotation may therefore 
be the result of close binary evolution (Pols et al. 1991). 

The visible companion (Section 2.2) seems to have 
colours consistent with a main sequence star making it 
unlikely to have been involved in the interaction that 
caused the rapid rotation. e Sgr is known to be an X-ray 
source (Berghoefer et al. 1996; Hubrig et al. 2001), which 
might indicate an unseen companion. 
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Figure 12. Modelled polarization for e Sgr at high spectral 
resolution. The pure Ohman effect (due to rotational line 
broadening) is seen most clearly on the Call K line at 393 
nm, the strongest metal line in this region. On the strongly 
Stark broadened hydrogen lines (the five strongest lines in 
this region) the Ohman effect in Q/I is swamped by other 
effects, but the Ohman effect is still present in U/I, and is a 
strong function of inclination 


5. FUTURE OBSERVATIONS 


Our results on e Sgr could be further tested by a num- 
ber of other observations. 


5.1. Interferometry 


Interferometry has been successful in imaging the ro- 
tational distortion of a number of stars of similar bright- 
ness and spectral type to € Sgr such as Regulus (Che 
et al. 2011), Altair (Monnier et al. 2007) and a Oph 
(Zhao et al. 2009). Such observations should be able 
to confirm the extreme rotational oblateness and de- 
tect the variation of surface brightness due to gravity 
darkening. Interferometry is better than polarimetry at 
measuring the inclination and hence distinguishing be- 
tween the range of models that are consistent with the 
polarimetry. 


5.2. Spectropolarimetry and the Ohman Effect 


Current observations of the polarization due to the 
rotational distortion of stars, such as those presented in 
this paper, have been made with broad-band polarime- 
ters because these instruments can achieve the high pre- 
cision needed to measure the small rotational signals. 
However, further information on stellar rotation could 


13 


be obtained from spectropolarimetric observations with 
sufficient resolution to resolve the lines. 

These effects include the Ohman effect first predicted 
in 1946 (Ohman 1946). The Ohman effect arises because 
at different points in the line profile of a rotationally 
broadened line the absorption arises in different parts of 
the star and therefore changes the integrated polariza- 
tion we observe. Models of the effect have been given 
by Collins & Cranmer (1991). The effect is small and 
has not yet been observed. Since e Sgr has the largest 
rotational signal yet seen in continuum polarization it 
is interesting to consider whether observations of the 
Ohman effect for this star are feasible and useful. 

Because our polarization modelling using SYN- 
SPEC/VLIDORT (Section 3.1) is carried out at high 
spectral resolution, our models predict what should be 
seen in spectropolarimetry. As in the modelling done by 
Collins & Cranmer (1991) the spectral lines are treated 
as pure absorbers. No line polarization processes are in- 
volved. The polarization is due to the same processes 
that produce the continuum polarization (primarily elec- 
tron scattering). 

Models for the polarization wavelength dependence of 
€ Sgr are shown in Figure 12 for the 380 — 450 nm region. 
This shows the problem with detecting the Ohman effect 
in a star of this spectral type. The spectrum is domi- 
nated by hydrogen lines, which are strongly Stark broad- 
ened by an amount that exceeds the rotational broad- 
ening. The pure Ohman effect is only seen on the metal 
lines, of which the Call K line at 393 nm is the strongest 
in this region. The large polarization features seen on 
the hydrogen lines in Q/I are not due to the Ohman 
effect, but are due to the fact that in these strong lines 
we are seeing different layers of the star in line and con- 
tinuum. 

The Ohman effect features seen in U/I are more in- 
teresting because they show a strong dependence on in- 
clination. The features seen on the hydrogen lines are 
smeared out by the Stark broadening, but may still be 
observable and would provide another way of refining 
the inclination of the star. Such an observation would 
require spectropolarimetry with R > 5000 and precision 
~ 5 ppm or better. 


5.3. Ultraviolet (UV) Polarimetry 


Models of rotating stars predict that polarization val- 
ues are expected to increase at shorter wavelengths 
(Collins et al. 1991; Jones et al. 2022) reaching ~1% 
(10000 ppm) at 120 nm for large w and high inclina- 
tion. Since the polarization due to a disk is not strongly 
wavelength dependent, and interstellar polarization de- 
creases at shorter wavelengths, observations in the UV 
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would remove any possible confusion between these dif- 
ferent sources. Such observations may be possible with 
future space missions similar to the proposed Polstar 
satellite (Scowen et al. 2022). 


6. CONCLUSIONS 


Multi-wavelength high-precision polarization observa- 
tions of e Sgr show the distinctive wavelength depen- 
dence expected for a rotationally distorted star, but with 
a much higher amplitude (~ 300 ppm) than seen in other 
rapidly rotating stars. In order to fit a rotating-star 
model, an additional positive polarization component 
is needed. We attribute this to polarization in a low- 
density (41071? g cm^?) equatorial gas disk which also 
produces the narrow shell absorption features seen in 
the spectrum. 

Detailed modelling of the polarization together with 
additional data including the spectral-energy distribu- 
tion and v, sin i reveals a very high rotation rate of w > 
0.995, the highest yet measured for a star in our galaxy. 
The modelling also results in a detailed set of stellar 
parameters. The star is rotating too rapidly to be con- 
sistent with single rotating star evolutionary tracks, sug- 
gesting that it is likely to be a product of binary evolu- 
tion. 

€ Sgr shares some features with Be stars in having a 
low density gas disk, weak Ha emission seen on one occa- 


sion, and multiple periodicities in its TESS light curve. 
However, the disk is of very low density and the TESS 
periodicites are of very low amplitude compared with 
typical Be stars. The results confirm that rapid rota- 
tion alone cannot be the cause of the Be phenomenon, 
and support the idea that an additional mechanism is 
needed to eject material into a Be disk. 

€ Sgr is an excellent target for future observations by 
interferometry, spectropolarimetry and UV polarimetry 
which would allow the extreme rotation to be confirmed 
and the stellar parameters refined. 
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APPENDIX 


A. INTERSTELLAR POLARIZATION 


In Figure 13 we present a map of the polarization of 
nearby stars. The stars have been selected such that 
their polarization is likely to be interstellar dominated. 
We interrogated observations made with high precision 
polarimeters, reported in Bailey et al. (2010); Cotton 
et al. (2016a); Marshall et al. (2016); Cotton et al. 
(2017a,b); Bailey et al. (2017); Bott et al. (2018); Cot- 
ton et al. (2019, 2020); Bailey et al. (2020a); Marshall 
et al. (2020); Piirola et al. (2020); Lewis et al. (2022); 
Marshall et al. (2023), and selected stars with spectral 
types between AO and K2.5 within 35°. From this list 
we removed rapidly rotating stars, those in close binary 
systems, magnetically active stars, systems hosting sig- 
nificant debris disks, and any other star known to display 
intrinsic polarization. To this we added a small number 
of interstellar controls newly observed for this work (Ta- 
ble 6). The observations were not all made in the same 
bandpass so, as we have done in past work, we assumed 
Amax = 470 nm and used the Serkowski Law to scale 


the polarization magnitude to correspond to 450 nm — 
representative of the g^ band. 

The right panel of Figure 13 illustrates that most 
stars near e Sgr have polarizations corresponding to < 
2 ppm/pc, though the closest star, HD 169586, is a lit- 
tle greater. This is in line with the trends previously 
reported in Cotton et al. (2017b), which if extrapolated 
predict 60.0 + 10.5 ppm for e Sgr (grey dot). 

The position angles of the interstellar controls are 
shown in the left panel of Figure 13. Cotton et al. 
(2017b) presents a formula for predicting position angle 
of interstellar polarization for a star based on the angu- 
lar separation weighted mean of nearby stars that gives 
163.1° + 29.1 (grey bar). However, this was only meant 
to apply to stars within 30 pc of the Sun, where po- 
tential controls are sparse. The paper also found better 
agreement than the nominal error between stars within 
10? of each other on the sky; the median position angle 
of the nearest 5 stars to e Sgr is 145.6?; it is likely that 
the interstellar position angle is within 15-20? of this. 


Interstellar Controls within 35.0 degrees of eps Sgr (HD 169022) 
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Figure 13. Plot of the polarization of stars near to e Sgr, representative of interstellar polarization, as projected onto the 
sky (left) — vectors represent the PA (clockwise N over E) — and debiased according to p = \/p? — p2 as a function of distance 
(right). Where multiple observations exist, the error-weighted average is displayed. The grey vector and dot represent the 
interstellar polarization of e Sgr as predicted by the model in Cotton et al. (2017b), whereas their black counterparts average the 
observations, and correct them for wavelength as if they are interstellar measurements. All observations are scaled to 450 nm 
as if they are purely interstellar with Amax=470 nm. The controls are numbered in order of angular separation as follows, 
1: HD 169586, 2: HD 165135, 3: HD 177474, 4: HD 166916, 6: HD 176687, 7: HD 165189, 8: HD 157919, 8 HD 156384, 9: 
HD 174309, 10: HD 160915, 11: HD 181240, 12: HD 160691, 13: HD 188114, 14: HD 151680, 15: HD 157172, 16: HD 191408, 17: 
HD 153580, 18: HD 192310, 19: HD 155125, 20: HD 151504, 21: HD 194640, 22: HD 180409, 23: HD 151192, 24: HD 165499, 25: 
HD 153631, 26: HD 146070, 27: HD 167425, 28: HD 197157, 29: HD 143114, 30: HD 173168, 31: HD 199288, 32: HD 164259, 33: 
HD 168723, 34: HD 162521, 35: HD 181391, 36: HD 197692, 37: HD 145518, 38: HD 192947, 39: HD 144766, 40: HD 185124, 
41: HD 177389, 42: HD 164651. Note that for stars with p < op the PA's in the left hand panels are not well defined, since 


Opa > 28°. 


Table 6. Observations of Interstellar Control Stars 


Control SpT Run UT Dwell Exp. Aeg Eff. q u p 0 
HD (9 (9 (m) (%) (ppm) (ppm) (ppm) (©) 
157919 F5II-IV 2018MAR 2018-03-27 19:22:38 1114 640 468.3 84.2 185+ 5.8 —20.4 + 5.6 27.5 + 5.7156.1+ 6.0 
160691 G3IV-V 2018MAR 2018-04-04 17:15:17 1780 1280 472.1 85.7 145+ 60 81+ 5.9 166+ 5.9 146+ 11.2 
165189 A6V 2018AUG 2018-09-01 10:51:27 971 640 464.4 58.7 —56.6 + 11.7 0.1 11.2 56.6 +11.5 899+ 5.9 
168723 KOIIT-IV 2018MAR 2018-03-29 19:00:51 1001 640 475.0 86.9 10.2+ 4.6 22.7+ 4.4 249+ 45 329+ 5.2 
2018JUL 2018-07-12 10:54:07 1100 640 475.3 86.1 17+ 46 170+ 45 1714 46 42.1 + 7.8 
177474 F8V 2017AUG 2017-08-14 13:49:16 1663 960 471.3 89.2 3.4 + 5.8 —17.6 + 5.9 179+ 5.8 140.5 + 9.9 
181240 ASIIT 2018MAR 2018-04-05 18:40:37 1716 1280 465.3 83.0 —88.4 + 6.9 —47.5 + 6.9 100.4 + 6.9 104.1 + 2.0 
2018MAR 2018-04-06 18:25:17 1726 1280 465.4 83.0 —104.2+ 6.7 —54.2 + 6.6 117.5 + 6.7 103.7 t 1.6 
191408 K2.5V — 2016JUN 2016-06-25 12:21:58 3029 1680 477.2 90.9  —8.6- 9.3 —32.0 + 8.7 33.14 9.01275+ 8.0 
192310 K24-V  2016JUN 2016-06-25 15:41:28 3464 2560 476.1 90.7 —9.8 + 7.2 —5.9+ 7.1 11.4+ 7.2105.5 + 22.3 
2018JUL 2018-07-20 13:21:57 2203 1600 475.3 86.0 —17.7 + 69 —3.5 + 8.3 180+ 7.6 95.6 + 13.9 
a 2018JUL 2018-07-23 11:07:49 1772 1280 627.0 84.4 —14.8 + 8.9—18.4 + 9.8 23.6 + 9.3 115.6 + 12.7 
a 2018JUL 2018-07-23 11:33:17 1051 640 627.0 84.4 14.0 Æ 12.0 —16.7 + 12.3 21.8 + 12.1 155.0 + 19.7 
192947 G9III 2018JUL 2018-07-18 12:25:25 1068 640 474.7 85.8 —16.1 + 4.7 15.7+ 4.7 2254+ 4.7 67.9+ 6.1 
Notes: The same aperture as used for the e Sgr observations specified in Table 1 in the same run was used. Spectral types 


(SpT) are from SIMBAD. P Indicates SDSS r’ filter paired with R PMT detectors; all other tabulated observations used the 


SDSS g' filter paired with the B PMT detectors. 
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